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We study a class of nonequilibrium lattice models on a ring where particles hop in a particular 
direction, from a site to one of its (say, right) nearest neighbours, with a rate that depends on the 
occupation of all the neighbouring sites within a range R. This finite range process (FRP) for R = 0 
reduces to the well known zero-range process (ZRP), giving rise to a factorized steady state (FSS) 
for any arbitrary hop rate. We show that, provided the hop rates satisfy a specific condition, the 
steady state of FRP can be written as a product of cluster-weight function of {R + 1) occupation 
variables. We show that, for a large class of cluster-weight functions, the cluster-factorized steady 
state admits a finite dimensional transfer-matrix formulation, which helps in calculating the spatial 
correlation functions and subsystem mass distributions exactly. We also discuss a criterion for which 
the FRP undergoes a condensation transition. 

PACS numbers: 05.40.-a ,64.60.De, 05.70.Fh 


I. INTRODUCTION 

Driven interacting many-particle systems [l[ have been 
of considerable interest in the past decades due to their 
rich transport properties, especially in lower dimensions. 
The zero range process (ZRP), a lattice gas model with¬ 
out any hardcore exclusions, is perhaps the simplest of 
them, which exhibit nontrivial static and dynamic prop¬ 
erties in the steady state. The ZRP was introduced 
as a mathematical model for interacting diffusing parti¬ 
cles and, since then, has found applications in different 
branches of science Si, such as in describing phase 
separation criterion in driven lattice gases network 
re-wiring [^, statics and dynamics of extended objects 
0 , [ 1 ], etc. Interestingly, the ZRP shows a condensation 
transition for some specific choices of particle hop rates 
for which, when the density becomes larger than a crit¬ 
ical density pc, a macroscopic number of particles accu¬ 
mulate on a single lattice site - representing a classical 
real-space analogue of the Bose-Einstein condensation. 
The ZRP has been generalized to multi-species models 
3, the misanthrc^ process [13] , urn models [HI ,the in¬ 
clusion process [I2| and inhomogeneous hop rates [l3| . 
etc. 

In the ZRP, the particles hop stochastically to one of 
the nearest neighbours with a rate that depends only on 
the number of particles on the departure site. As a con¬ 
sequence, the ZRP has a factorized steady state (FSS), 
which is amenable to exact analytic studies. However, 
when the hop rate depends on the neighbouring sites, the 
steady state does not factorize in general [ij, [T^l • In such 
situations, one may naturally expect a cluster-factorized 
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steady state (CESS), a straightforward generalization of 
the factorized steady state (FSS), where the steady state 
weight is a product of cluster-weight functions (see Eq. 
©) of several variables, i.e., the occupation numbers at 
two or more consecutive sites. 

In this paper, we study a class of nonequilibrium lattice 
models where particles hop in a particular direction, say 
from a site to its right nearest neighbour, where hop rates 
not only depend on the occupation of the departure site 
but also on the occupation of all of its neighbours within 
a range R; hereafter, we refer to this process as the finite 
range process (FRP). We demonstrate that, in one spa¬ 
tial dimension, one can have a CESS for various specific 
choices of hop rates; what we mean by the CESS here is 
that the steady state probability weight can be written 
as a product of functions of i?-|-1 variables, each of them 
being an occupation number in the cluster of R consecu¬ 
tive sites. A special case of the CESS with R = 1, called 
the pair-factorized steady state (PFSS), was recently pro¬ 
posed and studied in |14| where it was shown that, for a 
particular class of PFSS, the system can also undergo a 
condensation transition. Later, the PFSS has been found 
in continuous mass-transfer models [il[i3> in systems 
with open boundaries [l8l| and in random graphs 
etc. However, non-trivial spatial structure, which is not 
present in a FSS, has not been explored before. 

We show that, for a broad class of systems having 
a CESS with any R, there exists a finite dimensional 
transfer-matrix representation of the steady state. Being 
finite dimensional, these matrices are quite convenient to 
manipulate and help in exact calculations of spatial cor¬ 
relation functions of any order. Moreover, we propose a 
sufficient criterion for the hop rates that can give rise to 
condensation transition in FRP in general. Surprisingly, 
we find that a small perturbation to an FSS could destroy 
condensation transition, if any. 

The paper is organized as follows. In Sec. m we dis- 




2 


cuss the model and its steady state in general for the 
{R+ 1)- cluster. In Sec. Illli we formulate a transfer ma¬ 
trix method to calculate the correlation functions for 2- 
cluster and then, in Sec. IIVI we generalize the matrix for¬ 
mulation to i? > 1; continuous mass transfer models are 
also discussed in this section. Some useful applications of 
the FRP in the context of steady state thermodynamics 
for systems with short-ranged correlations, is studied in 
Sec. El Finally, Sec. EH provides conclusions, followed 
by open issues and discussions. 


II. MODEL 

The model (see Fig. 1) is defined on a one dimensional 
periodic lattice with sites labeled by i = 1,2,...L. Each 
site i has a non- negative integer variable rii representing 
the number of particles at that site (for a vacant site 
rii =0). Particle from any randomly chosen site i can hop 
to one of its nearest neighbours, say the right neighbour, 
with a rate that depends on the number of particles at 
all the sites which are within a range R with respect to 
the departure site: 

(. . . , Tli^ 7T,2-|_1, . . .) ^ (• ■ • ; ^2+1 “b 1, . • .) 

with rate u{ni-H ,..., rii,..., ni-t-fl)- (1) 

Clearly the total number of particles N = or the 

density p = N/L, is conserved by this dynamics. 


the factorized steady state can be achieved only for two 
cases - when the hop rate is u{k,m,n) = v{m) or when 
u{k,m,n) = w{m,n). The first case is the ZRP and the 
second one, where hop rate depends on both number of 
particles in both departure and arrival sites, is known 
as the misanthrope process (MP; see 0 for a review). 
Since a FRP with R > 1 includes i? = 1 as a special 
case, one expects that, except for the ZRP and the MP, 
there cannot be a factorized steady state (FSS) for these 
classes of systems. 

For the FRP, we first try whether a (i? -I- 1)— cluster- 
factorized form, 

L 

P{{ni}) oc ]^g( ni+R)5(^ni - N) (3) 

i—1 i 


with cluster-weight function g oi R + 1 occupation vari¬ 
ables, can be a steady state weight for Master equation 


|r((~.)) 


L 

u{ni-R, ..., rii,..., ni+R)P{{ni}) 

i=l 

L 


u{ni-R, ...,ni + l, rii+i - 1,.., Ui 
X P(... ,ni_i -b l,ni - 1,...). 


) 

( 4 ) 


Now, one can verify that a cluster-factorized form of 
steady state, as in Eq. ©, is indeed possible when the 
hop rate at site i satisfies the following condition 



i-R i-R+i i-i i i+i i+R-i i+R 


FIG. 1. (Color online) In the one dimensional finite range 
process (FRP) a particle hop from a site i to its right nearest 
neighbour with a rate that depends on occupation of site i 
(here m = 3) and all its neighbours within a range R. The 
lattice model, for a certain hop rate, can have an {R -\- 1)- 
cluster-factorized steady state. 

For R = 0, this model is identical to the zero range 
process (ZRP) Q with hop rate u{ni), an exactly solv¬ 
able non- equilibrium model that evolves to a factorized 
steady state (FSS) 

L 

P{{ni}) (xY[f{n^)6{Y^ni-N), (2) 

i—1 i 

with /(n) = 4"^® process considered in 

this paper is a generalized version of the ZRP and here¬ 
after we refer to it as finite range process (FRP). 

For i? > 0, the steady state of FRP in general cannot 
have an FSS as there are nonzero spatial correlations; 
however, there can be exceptions in specific cases. We 
provide explicit proof, in the Appendix, that, for i? = 1, 


R 

_ TT ff(hi—u+fc, hi-U-i-i-i-fc, . . . , hi-i-fe) 

k=0 diPi—R+k, Ri—R+l+k, ■ ■ ■ , ni+k) 

where fij = Uj — Sji. A simple way to prove this 
is to construct a pair-wise balance — for every hop 
that takes configuration C —>■ C' there is a suitable 
and unique configuration C” such that P{C)W{C —>■ 
C') = P{C'')W{C" C). For any configuration C = 

{..., ni_i. Hi,rii+i,..}, a particle hopping from site i can 
be balanced by taking C" = {..., rii-i-l-l, ni — \, ..} 

with hopping from i — I. Equation ([S]) is important as it 
says that any desired cluster-factorized state can be ob¬ 
tained in FRP by a choosing a suitable i?-range hop rate 
RiPi—R^ • ■ ■ ; • 1 Ri+r')- 

In the rest of the paper, we discuss various features of 
the cluster-factorized steady state and their applications. 

III. 2-CLUSTERS : PAIR FACTORIZED 
STEADY STATE (PFSS) 

Let us start with R = 1, for which the steady state 
is factorized as product of 2-site clusters, commonly 
known as the pair-factorized steady state (PFSS). In 
this case, particles hop from a site i to i -I- 1 with rate 
M(ni_i,rij,Tii+i) that depends on the occupation of de¬ 
parture site and its neighbours. To have a pair-factorized 
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steady state of the form 


P{{nz}) = ^ - Y[gin„ni+i)S(^ni-N) ( 6 ) 

-^L,N -7 


with a canonical partition function 

Zl,n = X! II ( X! ^ ) ’ 

{rii}i=l \ * / 

the hop rate must satisfy Eq. ([5]) with R= 1, 


For k = 1, we get density of the system p = (n) = 
InZL(z); by inverting this density-fugacity relation, 
one can express other observables as a function of p. 

This matrix formulation is quite general and works for 
any form of weight function g{ni, n^+i); however manag¬ 
ing infinite dimensional matrices is a challenging task. In 
the following, we show that, for a large class of weight 
functions, one can have a finite dimensional representa¬ 
tion which, in some cases, can even be extended to i? > 1 . 

Let us consider a weight function which has the follow¬ 
ing form 


, Tli , rij-l-l) 


g{n,_i,ni - 1 ) gjuj - 1 , rij+i) 

g{ni^i,ni) g[ni,ni+i) 


K 

-I-1) = ^ a„(ni)6K(nj+i), (11) 

K =0 


Unlike the FSS, the PFSS inherently generates spa¬ 
tial correlations and, like the FSS, it can lead to real- 
space condensation for certain hop rate H- This study 
has been lat er g eneralized on arbitrary graphs [T^ , open 
boundaries [l 8 | and for studying mass transport pro¬ 
cesses and condensation transition therein for discrete 
(particle) as well as continuous mass , etc. None of 
these studies, however, attempted to calculate the spa¬ 
tial correlations in these systems. In fact, the presence 
of spatial correlations can change the nature of transi¬ 
tions by creating spatially extended condensates with or 
without tunable shapes ( 20 |. 

To calculate spatial correlation functions we use the 
transfer matrix formulation which is possible for a large 
class of systems having a CFSS. For the purpose of illus¬ 
tration we mainly discuss this approach elaborately for 
the PFSS. Since the PFSS with any arbitrary cluster- 
weight function g{ni,ni+i) can be obtained from a suit¬ 
able hop rate u(ni_i,n^,n^+i) [as in Eq. ([71)], we rather 
focus on the functional form of g{ni,ni+i), not on the 
hop rate. In fact, any arbitrary function g{ni, n^+i) is an 
element of the infinite dimensional matrix 

00 00 

G 9 {n,n')\n){n'\ (8) 

n—0 n'—O 

where {jn)} are the standard infinite dimensional basis 
vectors which satisfy a completeness relation (njn') = 
Snn'- Then, in the grand canonical ensemble (GCE), 
where a fugacity z controls the density p, the partition 
sum can be written as 


ZUz) = Y, Zl,nz^ = Tr[T^] (9) 

N=0 


where the transfer matrix T has element (njTjn') = 
z("+" )/^p(n, n'). In the thermodynamic limit Zl{z) ~ 
^max (when Xmax, the largest eigenvalue of T is non¬ 
degenerate). Once we know the grand partition sum, we 
can calculate various observables; for example, all the 
moments for occupation number n at a site. 


(n^-) 


1 1 



Zl{z). 


( 10 ) 


where aK,(n),bi^{n) are arbitrary functions, not necessar¬ 
ily analytic. It is evident that g{ni,ni+i) can be written 
as an inner product of two {K + l)-dimensional vectors, 

g{ni,n^+i) = {a{ni)\P{ni+i)), ( 12 ) 

where 

{a{n )I = (ao(n), Oi (n),..., (n)) 

(^(n)| = ( 6 o(n), 61 (n),..., («))• (13) 

Then the partition sum in grand canonical ensemble is 
Zl{z) = Tr[T{z)^] with 

00 

r(z) = ^z"|/3(n))(a(n)| (14) 

n—0 

a {K + l)-dimensional matrix. Now the partition sum 
and the stationary correlation functions can be calculated 
easily. 

To illustrate this, let us consider a simple example by 
setting K = 1, 6 o(-) = 1 = ai(-)j renaming func¬ 
tions ao{.),bi{.) as /o(-)i/i(-) respectively. The weight 
function is now, 

g{ni,ni+i) = foiui) + fiim+i). (15) 


which we refer to as sum-form. This particular choice, 
i.e., a pair-factorized steady state with a weight function 
in sum-form, does not lead to condensation transition, 
which we discuss later in Sec. IV Bl Also, in Sec. IV Al we 
consider a general case of Eq. dm, which gives conden¬ 
sation transition, and we develop a possible criterion for 
the transition. 

For any functional form of /o(n) and fi{n) we always 
have an infinite dimensional representation given by Eq. 
([5]) . However, interestingly in this case, we can do away 
with the infinite dimensional representation and get a 
simple 2 -dimensional representation by taking, 

(a(n)l = (/o(n),l) and (^(n)| = ( 1 ,/i(n)). (16) 

The partition sum in GCE is then Z = Tr[T{z)^], where 


T{z) = Y^' 

n—0 


fo{n) 1 \ 

/o(n)/i(n) fiin) ) ' 


LZl{z) 


(17) 
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To see how the spatial correlation functions can be ob¬ 
tained, let us take a specific form of the functions /o(.) 
and /i(.), 


g{ni,ni+i) = 


+ 


{rii + ly {m+i + ly 


(18) 


where parameters v and 0 < g < 1 tune the hop rate 
of particles and g = 1 — g, corresponding to fo{n)/q = 
fi{n)/q = {n + In this case, the desired hop rate, 
for which the PFSS with weight-function as in Eq. ITSlis 
realized, is given by 


1, Tli , 



1" 

qn'^ + qirii-i + 1)'^ 

) 

q{ni + ly + g(ni_i -f ly 


q{n^+l + ly + qn'' 
q{ni+i + ly + qijii + 1)“' 


In the extreme limits g = 0 and g = 1, the model reduces 
to zero range process (details will be discussed in Sec. 

0. 

The transfer matrix , following Eq. CZl), becomes 


T(z )=1 {yy “fy,) 

z \ Li 2 v{z) quyz) ) 


(19) 


where Li^{z) are the Polylog functions. The eigenvalues 
of T are 


A± 


Liy{z) 

2z 



yq-qy 


4,qqzLi2u{z) 

(1 — z)Li^{zy 


( 20 ) 


The partition function Zl{z) = in the ther¬ 

modynamic limit {L oo) becomes Zl(z) ~ X^ and 
thus the density 


p(z) = z—\nX+. 
dz 


( 21 ) 


Throughout the paper, we calculate observables only in 
the thermodynamic limit. Let us consider q = \ and 
V = —1 (results for different g and v are discussed in Sec. 
0; here A± = ^(1 ± v^l + z)/{l — zY and the density is 


P = 


1 — z 2yi -I- 2 


( 22 ) 


Now we proceed to calculate the correlation functions, 
first the two-point correlation function and later the 
higher order. The two point correlation function is de- 


fined by 

C'(r) = (runi+r) - {ni){n^+r)- 

(23) 

For r > 0 

we have 



Tr[T'T^-^T'T^-^-y ^ 

y) - Tr[Ty ^ 

(24) 


where T' = dT/d{\nz). Eor q = \ and v = —1, we get 


C'(r) = 


ziy- 


4(l-f z)(l-z)2 




with ^ ^ = I In ^ I = I In | being the inverse cor¬ 

relation length. The correlation function for r = 0 is 
nothing but the variance cr^(p) of single-site occupation 
variable Ui, i.e.. 


CiO)^ayp) = {n^)-{ny^ 


Tr[T”T^-y 

Tr[r^] 


p2(26) 


where T" = d^T/d(]nzy and, again for g = 5, z/ = —1, 


C(0) 


z 

4(1 -2)2 


z2 -H 142-f 17 

(1 + ^) 


8 

\/l -f 2 


(27) 


Now, we turn our attention to higher order correlation 
functions. The 3-point correlation function, for example, 
is defined as 


C(ri,r2) — (^i) ) (^i-t-ri-t-r2 ) 


which, in terms of transfer matrix, can be evaluated from 
the expression 


C(ri,r2) =-- p .(28) 


We find that the three-point correlation function can be 
written in terms of the two-point correlation functions as 


C[T^,T2) = p [C{T^) + C{t2) - B{Z)C{T^)C{T2)] (29) 

where B{z) also depend on the parameters g and v; for 
g = 5 and i/ = —1, we get B{z) = 1-I-82(1-I-2)/(3-f 2)^. 
In a similar way, one can calculate all the higher order 
correlation functions exactly. 

To conclude, when the weight function g(ni,ni+i) is 
a sum of two functions as in Eq. (USD, the correlation 
length ^ = I In | remains finite at any density as A_ < 
A+ for any choice of g and v. 


IV. GENERALIZATIONS 
A. 3-clusters and general (R-(-l)-clusters 

In this section we consider some specific models of ERP 
with i? > 1 which give rise to {R + l)-cluster-factorized 
steady state. Corresponding partition function in the 
grand canonical ensemble would require contraction of a 
tensor product which is usually a hard task . Our aim 
here would be to obtain, if possible, a matrix formulation 
that can accommodate some cluster- factorized steady 
states for any R > 1. For R = 2 we have a 3-site cluster 
factorized steady state, 

= Y[g{ni,niy.l,n^^2)S n — . 

\i^l ) 

For illustration we consider a cluster-weight function, 

R 

R’i+l: • • ■ Ri+R^ — ^ ) f k^R’I+k) (^ 0 ) 

K — 0 


(25) 
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which is a simple generalization of the sum-form given in 
Eq. (fTSl) . We will now show that a grand partition func¬ 
tion of a finite range process which has a (i?-|- l)-cluster- 
factorized steady state with a weight function given by 
Eq. (1501) can be written as Zl{z) = Tr\T^] where 2 is 
the fugacity and T is a 2^-dimensional transfer matrix. 
Since we intend to obtain the transfer matrix for iter¬ 
atively, let us rewrite the transfer matrix given by Eq. 
(El) for i? = 1 in a convenient form, 


Ti(z) = ^z-Fi(n); F,{n) 

n—0 


fo{n) 1 \ 

/o(n)/i(n) /i(n) J 


(31) 

In a similar way, we extend to i? > 1 and write Tr = 
^"^=0 where the 2^- dimensional matrix can be 

written as 


^ { Fr_i Ar_iFr_i \ 

^ V JrFr-i /r^r-iFr_i J 


(32) 


2 

= [l - g{m^-2+k, m^+fe)”^] . (35) 

fc—0 

For small e, the model is equivalent to a discrete model 
where mass is measured in units of e. In fact, the resid¬ 
ual mass (actual mass modulo e) at any site does not 
change during evolution. The residual masses, each being 
smaller than a pre-defined value e which can be made ar¬ 
bitrary small, does not contribute to the asymptotic form 
of the hop rate. Thus we would obtain a transfer matrix 
T 2 discussed in the previous section, with /o,i, 2 (™) = fn, 
but the sum will now be replaced by an integral 
J dm. Defining, a chemical potential fj, (where 2 : = e^), 
we get the transfer matrix, as in Eq. (1331) . 


F(m) 



Ai 0 0 \ 

1 1 fj. 

1 0 0 




1/ 


(36) 


using a constant matrix 


Ar 


0 0 \ 
/ 2 H -1 0 J ’ 


where / 2 H -1 is a 2^“^-dimensional identity matrix. For 
i? = 0, we take Aq = \. Since R = 0 corresponds to 
the ZRP which has a factorized steady state, we have 
Fo(n) = fo{n), which is a scalar. Clearly Fi in Eq. (1^ 
satisfy Eq. (1321) . A little more algebra would show that 
the transfer matrix for i? = 2 is 


This matrix has eigenvalues ^{A, Aie^*®, A 2 }, where A 
(the largest eigenvalue), Ai, 6 and A 2 are independent 
of g, and their approximate numerical values are A ~ 
3.86841, Ai ft! 1.10465, 9 a; 1.87254 and A 2 ~ -0.21184. 
In the thermodynamic limit, the partition function is 
Z]^ = {\/, and density p = —2/^4. The two-point 
correlation function for r > 0 is 


C(r) = {mimi+r) - 
A 2 


= P 


C2 


A 


■ 2ci 


cos(r0 -I- a) 


(37) 


n 


F I T^iF I 

/ 2 F 1 /2^iFi 


^2"F2(n). (33) 


From the transfer matrix, one can, in principle, calcu¬ 
late the expectation value of any desired observable. We 
will not discuss further the finite range process i? > 1; 
the finite dimensional transfer matrix is expected to gen¬ 
erate spatial correlations which was absent in the ZRP. 
We discuss some of the models in details which undergo 
condensation transitions (see Sec. 0. 


B. Continuous mass model 

Until now, we have studied CFSS on a one dimensional 
lattice with each site having a discrete variable, called the 
occupation variables or number of particles. The model 
and the matrix formulation can be extended, without any 
particular difficulty, to systems with continuous mass m. 
As an example, let us consider 

5(TOi,mj+i,mj+2) =mi + mi+i +mi+2- (34) 

A 3-cluster-factorized steady state with the above weight- 
function can be obtained when e amount of mass is trans¬ 
ferred from site i to z -I- 1 with the rate 


where, ci = 0.3380, C 2 = —0.0375 and a = 0.1804. And, 
for r = 0, the correlation (actually cr^) is <7(0) = (m^) — 
= 0.6704p2. 


V. APPLICATIONS 


A. Condensation transition 


One important feature in these simple one dimensional 
models is that they can exhibit condensation transition 
at a finite density when one or more parameters in the 
rate functions are tuned. To demonstrate the possibility 
of a condensation transition in the CFSS, for any R, we 
consider the weight of (R + l)-cluster to be. 


^(77.2, 7 T,j_|_l, ..., ?T. 2 _|_r) — 


^j=o ''^7+i 


1 K 




(38) 


where K, v and q are positive and K is an integer. This 
steady state weight can be generated from a hop rate 
given by Eq. (ED, 



1 

q + IT'i-j+k 



u{mi-2, Wj- 1 , mj, Wj+i, TOi+2) 


( 39 ) 
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Case with R = 1 (PFSS). — We first consider R = 1. 
It is easy to see that for any K, the weight function Eq. 
([Ml) can be expressed as Eq. (HU with suitable choice of 
anin) and where k varies from 0 to K, leading to 

a {K + 1) dimensional transition matrix. We further set 
the parameters K = 1 = q\ this gives rise to a PESS, as 
in Eq. (jH), with g{ni ,, rii+i) = {m + m+i + l)/(ni + 1)", 
which can be realized when a particle hops out from a 
site i (to the right neighbour), having rii > 0 particles, 
with the following rate 

iV n^ + m+i 

rij / l + ni + ni-il + m + rii+i 



For this case, we can obtain exact results follow¬ 
ing the matrix formulation developed here. Eirst 
we write g{m,n) = {a{m)\j3{n)) where {a{m)\ = 

((to-I-I)”"", (m-I-1)^““^), (/3(n)| = (n, 1). Thus the 
grand partition function can be written as Z(z) = 
Tr{T^) with 


T = ^|/3(n))(a(n)|z" 

n—0 


1 / Liu-i{z) Liv-2{z) \ 

2 Li^{z) Li^_i{z) j 


The eigenvalues of T are 

A±(-2) = ^ (^Tv_i(z) ± Li„{z)Li^_2{z)^ , 

which leads to the density-fugacity relation p{z) = 
z\'Az)/\+{z) and the critical density pc = lim p{z). It 

turns out that for < 4, pc diverges — indicating a flnid 
phase for any density. For i/ > 4 we get, 

26(0 + 2C(i^ -1)06(0 06(0 + C(i^ - 1) 

where 6(*^) = Ci^)Ci^ ~ s-nd Ci^) are Riemann zeta 
functions. Thus, for z/ > 4 we have a condensate when 
density exceeds this critical value. Unlike the ZRP, where 
particles at different sites are not correlated, here we have 
non-vanishing correlation that extends up to a length 
scale 6-2^) = I la h(z) l~^ which is finite throughout. 

Case with R> 2 (CFSS). — It is straightforward to ex¬ 
tend the matrix formalism to i? > 1 when K = 1. First, 
let us take = 0. In this case, the weight function g takes 
a sum-form given by Eq. O, for which we have already 
constructed a general transfer-matrix. For z/ > 0, the 
dimension of the transfer matrix remains the same as in 
zz = 0; it is only that each element of F/j in Eq. (IMl) will 
be multiplied by an extra factor (rii -|- 1)“^. We omit the 
exact analytic expressions of the density-fugacity relation 
and the critical density - the calculations are straightfor¬ 
ward but the expressions are very long. Only the numer¬ 
ical values of critical density are tabulated in Table U for 
different parameters. 

Criterion for condensation transition.— For the ZRP, 
it is well known that, provided the hop rate uo{n) has an 
asymptotic form 


uo{n) = I -I- — -I- ... 
n'^ 


(41) 


TABLE I. Critical density pc for K — 1. 


p 

9 = 1 

R = 1 

9 = 1 

R = 2 

9 = 2 

R=1 

5 

0.3254 

00 

0.1591 

6 

0.1054 

0.2773 

0.0544 

7 

0.0429 

0.0981 

0.0228 


condensation occurs at a finite density, when cr < 1, or 
when tr = 1 but b > 2. This criterion can be extended to 
any other system (without any constraint on occupation 
number) when the steady state has a factorized form ([2]); 
one needs to consider and effective rate function uo(n) = 
f{n—l)/f(n) and find its asymptotic form. This criterion 
determines whether a model can undergo a condensation 
transition and helps in understanding phase coexistence 
in hardcore lattice gas models [1, . 

Such a criterion for cluster-factorized steady state 
would be very useful. At present, we do not have a gen¬ 
eral criterion, but the examples studied above suggest a 
sufficient condition for CFSS to have condensation. If 
the rate function can be expanded as 

/ \ \ ' d3jy(^ni—ni—i^ ^ 

.., = y - , 

n K 

i /=0 ‘ 

the condensation transition can occur when both the con¬ 
ditions 


(i) both Bq and are constant 
(R)Bi/Bo > 2 (42) 


are satisfied. This is only a simple generalization of the 
criterion of condensation in the ZRP. Effectively, Bi/Bq 
plays the role of b in Eq. 62). As the hop rate in Eq. 
can be expanded as 


u {... ,ni_irzi,nj+i ...) = 1 -|- 


i2-K{R+1) 


oi-), 

nf 


and thus Bq = 1 and B 2 = i'—K{R+l), the criterion cor¬ 
rectly predicts the condensation which occurs only when 
zz > Ar(i?-|-l)-|-2. This is same as the usual condensation 
criterion in the ZRP if we treat b = B 2 / Bq. In this par¬ 
ticular case, we have also checked that moments (n^) as 
a function of z, in leading order, are the same as that in 
the ZRP with corresponding b (see Eq. (HTl) ). This crite¬ 
rion, however, cannot be applied to some of the following 
cases studied recently, such as, the misanthrope process 
[iU| and the PFSS [ij . For the first case, Bq and Bi are 
not constants and, for the later case, hop rates are not 
analytic functions. A criterion of condensation, which 
can apply to a cluster-factorized steady state in general 
is desirable and remains a challenge. 

We end this section with the following remark. The 
condensation transition here is different from that ob¬ 
tained for PEES by Evans et. al. [l^. There, one ob¬ 
serves an extended condensate where both the size and 
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FIG. 2. (Color online) Particle distribution in FRP with 
weight function (I38II after t = 10® MCS, starting from a ran¬ 
dom distribution of particles. Density is p = pc + 0.01. The 
critical density for 7Z = 0 is pc = 0.01925; the same for R> 0 
are taken from table |T] Clearly, for all cases, the condensate 
is localized to a single site. The condensate size is written 
beside the condensate site. 



the spatial extent of condensate scales with system size 
as '/L. This indicates that the transition is associated 
with a diverging spatial correlation length. Whereas for 
the PFSS (and the CFSS) studied here, the correlation 
length remains finite throughout and the transition is 
characterized by a diverging mass fluctuation, as in the 
ZRP. The condensate is also localized to a single site (see 
Fig. [5]). A detailed comparison of nature of condensate 
would be reported elsewhere (^ . 


sum-form. To illustrate this, we consider a simple zero 
range process with weight function f{n) = l/(n -I- 1)'^ or 
hop rate u{n) = f{n — l)//(n) = {n + ly /n" and add a 
perturbative term get a new weight function 

gijii, m+i) = (1 - q)f{ni) -b qj{n^+i) (43) 


which depends on occupation of two consecutive sites. 
Here Q < q < l,q = 1 — q and we choose f{n) = l/(n -b 
1)“^. A pair-factorized state, as in Eq. ([5]), with the above 
weight function occurs when particles hop rate is 




_ g/(n»-i) + g/(n» - 1) 
qf{n^-i)+qf(ni) 

^, qf{ni - 1) + qf{n,+i) 
qf{n^) +qf{n^+i) 


For both q = 0 and g = 1 we have a factorized steady 
state, as in Eq. which corresponds to the ZRP with 
particle hop rate 


u{n) 


f{n-y 

f{n) 



:^l + --bO(4)- (44) 

n 


Thus we expect a condensation transition for <7 = 0,1 
when u > 2 and the density is larger than a critical value 
Pc- In this case the Z{z) = F{zy (the transfer matrix 
T(z) reduces to a scalar), where F{z) = J2‘^=o h{n)z'^ = 
Lii,{z). The density is p = z-^F{z) and thus the critical 
density for g = 0,1 is 


Pc = lim p{z) 

z —>1 


oo for v < 2, 

^fer-1 forJ^>2. 


(45) 


B. Pair factorized state with weight function in 
sum-form 


In this section, we first show that a pair factorized 
steady state with weight function < 7 ( 71 ,, rii+i) = f^ini) -b 
/i(ni+i), which we refer to as sum-form, cannot give rise 
to condensation. Then, we demonstrate this considering 
a perturbation to a ZRP that converts the existing fac¬ 
torized steady state of the ZRP to a PFSS with weight 
function in the sum-form. For the PFSS with weight 
function in the sum-form, the transfer matrix T(z) is 
given by Eq. (HZl). 

The largest eigenvalue of the matrix A+ = ^(Tfi -b 
T 22 + y{Til + 7 ) 22 )^ — 427, where V is the determinant 
of T can be used in Eq. (ED to get the density p{z). 
With some straightforward algebraic manipulations, one 
can show that the maximum density at z = Zc = ^ is, 


Pc = 


lim p{z) = liin - 

Z^l 2—>-1 Z 


1 dT2i 1 
T 21 dz 1 — z 


Clearly pc diverges independent of the first term, leading 
to a conclusion that there can not be a condensation 
transition at any finite density. Thus, a PFSS cannot 
have condensation transition if the weight function has a 



z 


FIG. 3. (Color online) Small perturbation to the ZRP: For 
small q = 10“^, 10“® or 10“^, density p{z) diverges when 
z 1. However for g = 0 or for g = 1, pc = p(l) = ~ 1 

is finite, leading to a condensation transition when p > pc. 
Inset shows the phase diagram for 1 / = 3. 

The phase-diagram of the condensation transition in 
the p-v plane is shown in Fig. [31 The critical line Pc 
separates the condensate phase from the fluid phase. For 




































a general 0 < g < 1 , we need to calculate the density 
using Eqs. (HU), (EOl), and (EB), 

p{z) = - f. \r - / - zfLi^{z)- 

a[a — (1 — z)Lii,[z)) 

(1 - z)a) + 2 qqz{Li 2 u{z) + (1 - z)Li 2 u-i{z))] — 1 

where _ 

a(g, z) = ^J{q- g)2(l - zYLil{z) + ^qqz{l - z)Li 2 v[z). 


In the limit z —>■ 1, p{z) diverges for all z/ > 0 and thus 
the condensation transition is destroyed. It is somewhat 
surprising why for any non-zero q however small, the con¬ 
densation transition is destroyed. It seems that this per¬ 
turbation, which takes the factorized steady state of the 
ZRP to a pair-factorized steady state, forces the conden¬ 
sation to disappear. One could understand this following 
the criterion (|42l) . For zz = 3, the rate for general q has 
an asymptotic form (be. when rn —> oo) ) 


= 1-1-3 


g^(l -b -b g^(l -b m+iY 

qqnj 


Thus, here Bi/Bq = 0 and therefore we should not expect 
condensation for this hop rate. It can be shown easily 
that for any zz > 1 the asymptotic form of the hop rate 
does not satisfy condition {ii) of ansatz ruling out 
the possibility of a condensation transition. 


C. Subsystem mass distribution 

It was argued in recent works [H, [1^ [2^ that, for 
systems with short-ranged interaction, irrespective of 
whether they are in equilibrium or not, one could ob¬ 
tain a state function which plays the role of a free en¬ 
ergy function. It was shown in [2^ that the steady 
state distribution (m) of mass in subsystems of volume 
V ^ ^ can be determined from the functional dependence 
of the scaled variance cr'^{p) = {{rn?) — {m)'^)/v, in the 
limit of large v, on the mass density p. When a^{p) oc p^ 
is a quadratic function of density p, the subsystem mass 
distribution can be characterized through a gamma dis¬ 
tribution, i.e., Pv{m) oc exp(pTO), where p = —rj/p 
is an equilibrium-like chemical potential. The exponent 
rj however depends on the details of the model and it can 
be calculated from the knowledge of two point-correlation 
function only. The matrix formulation developed here for 
the CFSS can thus help in determination of p. 

To illustrate this, let us consider a continuous hnite 
range process with i? = I, and calculate explicitly the 
variance of the subsystem mass. Consider the following 
homogeneous weight function for a pair-factorized steady 
state. 


g{m^, m^+i) = ml + c m]ml_^_^ (46) 

The grand partition sum is Z = rr[T^] where the trans¬ 
fer matrix T(p) (p = ln(z) is the corresponding chemical 


potential) is given below 


Tip) 


1 


r (5 + i) \ 

cr(5 + 1 )) 


(47) 


where r(.) are Gamma functions. Eigenvalues of T(p) 
are A± = A±(A, 7 , c)/p^+^ where 

2 A± (<5, 7 , c ) = (1 -b c)r((5 -b 1 ) _ 

± + 1)^(1 - c)2 -b 4cr(2(5 - 7 -b l)r (7 -b I) 

and the particle density is 


P = 


dp 


lnA+ = — 


^+1 

P 


(48) 


implying a fluctuation-response (FR) relation 


I = "'(z). («) 

analogous to the fluctuation-dissipation theorem in equi¬ 
librium. Now, as shown below, one can check the above 
FR relation by explicitly calculating both sides of Eq. 
(|4^ . The r.h.s of Eq. (|4^ can be calculated by integrat¬ 
ing two-point correlation function cr^(p) = 
using Eq. (|24)) . 

C{r) = {rurii+r) - p^ = p^ [A{r) - 1] (50) 


where, for r > 0 , 

A{r) = I -b 



(^- 7 ) 2 /(^- 1)2 

_ r(j+i)^ 

r(25-7-i-i)r(7-i-i) 


and 

^ ^ 2c (A- 7 )^r( 2 A -7 + I)r( 7 +l) 

^ ’ 5+1 ^A+(<5 + 1)2 2A+-(l + c)r(A + l)- 


In this system, the gap (A+ — A_) between the two 
eigenvalues is nonzero and the correlation length ^ = 
|ln^|“^ is finite. Therefore, following Ref. [l^, the 
subsystem mass distribution Py{m), for c ^ is a 
gamma distribution where the exponent p can be written, 
using Eq. m , as 


r——oc) 


Note that the exponent p depends only on the homogene¬ 
ity exponent 6 but neither on 7 nor on c. The left-hand 
side, the compressibility dp/dp, of Eq. (H51) gives the 
same p = P^{^)~^ = 5 + 1, by differentiating the ex¬ 
pression p = —(5 + l)/p in Eq. (H51) with respect to p; 
this is a proof that the fluctuation-response relation in¬ 
deed holds here and also is consistent with the additivity 
property proposed earlier for these systems [H, [l^l . 

In principle, the single-site mass distribution (for v = 
1 ) can be calculated straightforwardly from the moments, 
but the exact closed form expression is hard to obtain. 
In this regard, this formulation [ 2 ^ [2^ for obtaining the 
subsystem mass distribution from the two-point correla¬ 
tion function is quite useful in obtaining the macroscopic 
behaviour of the system. 
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VI. SUMMARY 


APPENDIX 


We have introduced a class of nonequilibrium finite 
range processes (FRP) where particles on a one dimen¬ 
sional periodic lattice can hop in a particular direction, 
from one site to one of its nearest neighbours, with a rate 
that depends on the occupation of all the sites within 
a range R starting from the departure site. We show 
that, for certain specific functional forms of the hop rates, 
the FRP has a cluster-factorized steady state (CFSS), 
i.e., the steady state probability of a microstate can be 
written as a product of cluster-weight functions g having 
{R+ 1) arguments - the occupation numbers of {R+ 1) 
consecutive sites. The model with i? = 0 reduces to the 
familiar zero range process (ZRP), which has factorized 
steady state. 

The CFSS with R = \ reduces to the pair-factorized 
steady state (PFSS) and its steady state can always be 
represented by an infinite-dimensional transfer matrix. 
However, for the CFSS with i? > 1, a matrix formu¬ 
lation is not guaranteed. In this work, we show that, 
for a large class of systems having CFSS with i? > 0, 
there exists a finite dimensional matrix representation. 
Being finite dimensional, these matrices are easy to ma¬ 
nipulate and thus help in exactly calculating the n-point 
correlation functions for any n. The two-point corre¬ 
lation function (n = 2) can be utilized to characterize 
the subsystem mass distribution in these nonequilibrium 
systems in terms of a nonequilibrium chemical potential 
and a free-energy function, which are obtained through 
a fluctuation-response relation [l^, - analogous to the 

equilibrium fluctuation-dissipation theorem. 

Even though the transfer-matrix is finite dimensional, 
the CFSS can undergo a condensation phase transition. 
We obtain a sufficient condition for the condensation 
transition for a particular class of hop rates in the FRP 
in general. The nature of the condensation transition 
studied in this paper are however different from those 
studied in systems having a PFSS [l^; the condensate 
here remains localized, as in the ZRP, in contrast to the 
extended condensate observed in Moreover, 

the condensation transition studied here occurs solely due 
to the diverging mass fluctuations at certain critical den¬ 
sity Pc, not due to a diverging correlation length; in fact, 
the correlation length remains finite throughout. 

We should mention that it is always possible to con¬ 
struct a hopping dynamics of FRP so that it evolves to a 
desired steady state which is cluster-factorized. However, 
for a given hop rate, there is no simple way to check if 
it can give rise to a CFSS. In our opinion, the FRP is a 
very general class of models as it includes the Ising model, 
Potts model, misanthrope process, urn models, symmet¬ 
ric and asymmetric exclusion processes on a ring and 
many other models (one of them, of course, zero range 
process). More importantly, the method developed here, 
could help in finding the exact steady state structure in 
models even when the interactions extend beyond two 
sites. 


In this Appendix, we provide an argument that FRP 
can have a factorized steady state only for i? = 0 (namely 
the ZRP) and for some specific misanthrope process (spe¬ 
cial cases oi R= 1). For any i? > 1, however, one cannot 
have a factorized steady state in general. First we con¬ 
sider i? = 1 and show that, in this case, the hop rate 
reduces to those in the ZRP or in the misanthrope pro¬ 
cess, when one demands a factorized steady state. One 
can construct a general proof in a similar way, that con¬ 
dition of FSS would reduce the hop rate of FRP with 
i? > 1 to the ZRP or the misanthrope process. In Ap¬ 
pendix B, we provide a proof of the above for the hop 
rates which can be written in a product form. 


APPENDIX A : NO FSS FOR R = 1 


In this section, we show that, for i? = I, one cannot 
have a factorized steady state for the general hop rate 
M(ni_i,ni,ni_|_i). The Master equation for FRP for gen¬ 
eral i? > 0 is 




ni,...,ni+R), (52) 


where 


F{ni_R ,.., m ,.., Ui+n) = u{ni-R ,.., ni ,.., ni+R)P{{ni}) 

uijli—R^ .., Tli -t- I, 1, .., 

X P{..,ni + l,ni+i-l,..). (53) 

In the steady state, right hand side of Eq. (15^ must 
vanish, which can happen if 


R iFi— Ri * ■ * , , • ■ • , Ri+Ft) — hijli— R^ • ■ • , ^2, • ■ • , Fii+R—l) 

- h{ni^R+i,... .m,... ,ni+R) ( 54 ) 


for some arbitrary function h of 2R arguments. Note, 
that the above cancellation scheme is only a sufficient 
condition. 

Now let us consider i? = 1, and demand that the steady 
state has a factorized form given by Eq. (I2|). Then 


u{ni-i,ni + 1, Ui+i - I) 


f{rii + 1 ) firii+i - 1 ) 


f{ni) f{n^+l) 

-u{ni-i,ni,ni+i) = h{ni,ni+i) - h(ni_i,ni) (55) 


where h is an arbitrary function, yet to be determined. 
Since the hop rate u(ni_i,n^,n^+i) = 0 when Ui = 0 
and we must have a boundary condition f{m < 0) = 
0, we can use specific values of n^s in Eq. (1551) to find 
recursion relation for h. For = 0 = rii+i equation (l55l) 
in h{ni-i,0) = /i(0, 0). Again putting rii+i = 0 = Ui-i 
we get h{0, m) — h{0, 0) = ^(0, m, 0). 

These two conditions leaves Eq. ([55]) for rii = 0 as 

/ 1 /K-ti - 1) rrs n\ 

/(n-+i) ^ m( 0, "raz-n, 0)- 
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Clearly, in order to be consistent, 1,rii+i) must 

be independent of rii-i. For convenience, without any 
loss of generality, lets set u(ni_i, l,ni+i) = u(0, l,ni+i). 
Thus, to have the factorized steady state for i? = 1, the 
hop rate rii+i) must satisfy 


u{ni-i,ni + l,nj+i - 1) 


M(0,l,ni) 'u(0,ni+i,0) 

u(0, n^ + 1,0) m(0, 1, Tii+i - 1) 
= u(ni,ni+i,0) - it(ni_i,ni,0). 


Now if we take Ui = 1 and use u(ni_i, 1,rii+i) = 
m( 0, 1, rii+i) in the above equation to rearrange the terms, 
we have 


u{ni-i,2,ni+i - 1) 
-m( 0, l,nj+i) = 


u(0,l,l) u{0,ni+i,0) 

u(0, 2, 0) u(0, l,ni+i - 1) 
= u(l,nj+i,0) - u(0,1,0), 


which implies that u{ni-i,2,ni+i) must be indepen¬ 
dent of rii-i. A similar recursion would result that 
u{ni-i, rii, rii+i) must be independent of n^-i. This again 
reflects the fact that a factorized steady state is possible 
for R = 1 only when hop rate is m = u{ni,ni^i) i.e. the 
process is a misanthrope process. 


APPENDIX B : NO FSS FOR HOP RATE 
HAVING PRODUCT FORM 

In this section, we show that the FRP, for any i? > 0, 
cannot have a FSS when the hop rate has the following 
product form. 


— ..., Ui+Bi) 


(57) 


where Vk = Vkijii+k) and 


G(ni,nj+i) 


- uo(nj -I- l)vi{ni+i) 

+ vo{ni)vi{ni+i). 


fjnj + l)f{ni+i - 1) 


Now differentiating both sides of Eq. (l57l) with respect 
to Ui-H and Ui+R, we have 


---- V-R+i. ..V-1V2 ■. .VR-iG{ni,ni+i) = 0 . 

OUi-R OUi+R 

This implies that, either V-Rijii-R) or VR(ni+R) must be 
a constant, because the other solution /(n) = l/voin) = 
l/vi{n) cannot be accepted as it means vi(0) = r'o(O) = 
0, ie., a particle cannot be transferred to a vacant neigh¬ 
bouring site. So, let us proceed with v-r =constant (say 
1). Then Eq. (15^ reads as 

vi-R..V-iV2:VRG{ni, ni+i) = h{ni_R, .., Ui, .., m+R-i) 
h[ni—R-{-i^ .., Tii^ .., riiJf-R'). 

Clearly for this equation to be valid its right hand side 
must be independent of ui-r and that in turn leads to 
h{xi,X2, ...,Xk)= h(x2, . . . ,Xfc). 

This way we can eliminate one variable at each step 
until finally we reach to 


R 

u{n,_R,...,n^,...,n^+R)= Vj{n,+j). (56) 

3=-R 

The Master equation along with a demand of a factorized 
steady state of the form and then Eqs. (l551l and (IMl) 
together, implies 

V-R ... V-1V2 ■ ■. VRG{n^, Ui+i) = h{n^-R ,..., n^+R-i) 


v{ni-R, ..., rii, ..., ui+r) = v{ni, n^+i) = vo{ni)vi{ni+i), 
j , IN / + 

and vo{ni -I- l)'Ui(ni+i - 1)———7 -—-— 

JKn,) f{ni+i) 

-vo(ni)v-i(ni+i) = h(ni+i) - h{ni). 


This is exactly the criterion for having a factorized steady 
state in misanthrope process with a hop rate that has a 
product form u{ni,ni+i) = vo{ni)vi{ni+i) 0- 
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